Entanglement magnification induced by local manipulations 
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We study the entanglement capability of the evolution of a pair of qubits evolving under unitary 
dynamics, when the local dynamical parameters cannot be modified during the time-evolution. 
Unlike the fast local control regime, we find that local and non-local contributions to the dynamics 
are strictly interconnected. Moreover, it is possible to strongly increase the entanglement capability 
by suitably initializing the characteristic energies of the two parties. 
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Introduction. — A basic ingredient for the implemen- 
tation of quantum technologies is the ability to control 
some fundamental processes involving a pair of qubits, 
that is, the capacity of influencing the dynamics of this 
system through external actions. When the interaction 
with the environment is negligible, the system dynamics 
is given by the family of unitary transformations Ut gen- 
erated by the total Hamiltonian Ht = Hi + II2 + Hj, 
where Hi, H2 are one-qubit contributions, and Hj rep- 
resents the interaction. This term is responsible for the 
generation of the peculiar quantum correlations called 
entanglement. The growing interest in the quantum the- 
ories of information and computation is largely due to 
these correlations, and to their potential applications. 
Therefore, it is of fundamental relevance to character- 
ize the entanglement, to understand its properties and 
to study the processes leading to its creation. The rel- 
evant quantities involved are determined by the specific 
physical realization of the qubits pair. 

The standard control setting is the so-called Local Uni- 
tary control (LU), where the control operations act lo- 
cally (they influence Hi and H2, not Hj). Moreover, 
motivated by the available technologies, it is usually 
assumed that these operations are performed instanta- 
neously P, 01 (e.g., by means of short laser pulses). An 
arbitrary unitary transformation Ut can be written as 



Ut = An(E> Bne 



-iHt„ 



-iHt-2 



Ai ® Bie 



Aq ® Bo, 
(1) 

representing a succession of entangling evolutions, acting 
for time intervals . . . i„, interspersed by instantaneous 
local evolutions, represented by the A^ and Bi operators, 
that can be arbitrarily manipulated. In general H ^ Hi, 
since the physical interaction can be used to simulate the 
dynamics generated by a different non-local contribution 
[21 • In (P), the non-local and local parts are mutually 
independent. 

Other control schemes have been developed for unitary 
dynamics, inspired by different experimental scenarios. 
In particular, the results presented in this work are rel- 
evant for the indirect control methods, in which an aux- 
iliary system (ancilla) is used to manipulate the target 
system through their mutual interaction 0] ■ 



In the framework of LU control, the processes involv- 
ing the qubits pair have been widely investigated in the 
past years. For fast local controls, the time-optimal gen- 
eration of entanglement has been considered in us- 
ing the idea of entanglement capability of the interaction. 
For an arbitrary unitary operator, the maximal achiev- 
able entanglement, and the corresponding uncorrelated 
initial state, have been characterized in [6|. A geomet- 
rical characterization of the entangling gates has been 
given in in terms of the coefficients describing the 
entanglement production. Other relevant topics, not di- 
rectly connected to the entanglement generation, have 
been discussed in the literature, as the simulation of non- 
local gates [H, d, H, S Il3| 7 and the ability of transmitting 
classical as well as quantum information Q. 

Even if usually well justified, the expression ([1]) is only 
an approximation, since it involves instantaneous actions. 
In this work, we relax this assumption by considering the 
opposite regime, in which the controls cannot be modified 
during the time evolution. In this case, it is not possible 
to distinguish between independent local and non-local 
contributions, as in ([T]). It turns out that the entangling 
part of the dynamics inherits an explicit dependence on 
Hi and H2 as long as [Hi + H2,Hj] ^ 0. Local control 
and entangling dynamics become deeply interconnected, 
and this dependence can be used to manipulate the non- 
local part by means of local control. In this letter we 
derive the entanglement capability in this regime, and 
show that particular choices of the local parameters can 
highly increase it. 

Our analysis complements previous results obtained 
under the fast LU control assumption. Moreover, it is of 
interest in the context of indirect control methods, where 
the entanglement between target and auxiliary systems is 
fundamental to perform manipulations, and the param- 
eters of the ancilla have to be fixed accordingly. 

Cartan decomposition of unitary operators. — We as- 
sume that the time-evolution of the pair of qubits is given 
by a family of unitary maps generated by the total Hamil- 
tonian. A fundamental property, known as Cartan de- 
composition (or canonical decomposition), factorizes an 
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arbitrary unitary operator (7 S i7(4) as 



U = LAK, 



(2) 



where L,K ^ SU{2) SU{2) are local operators, and 
A = e"", with a an element of the Cartan subalgebra of 
the Lie algebra u(4) = span{iaj (g) aj,j = x,y,z), CFj 
the Pauli operators for the two subsystems. The only 
operator that can correlate the two systems is A, that is 
called the entangling part. L is irrelevant when dealing 
with entanglement generation. 

This decomposition clarifies the role of the 16 param- 
eters entering an arbitrary unitary transformation J7(4): 
the local contributions are characterized by twelve of 
them, the composite evolution by three of them, and fi- 
nally there is a not-relevant overall phase. In particular. 



(3) 



where Oj {j = x,y,z) and \j {j = 1,...4) are real 
constants embodying the entanglement capability of the 
channel. In Q, we have introduced the eigenvalues \j 
and eigenvectors \j) of a (the so-called magic basis 
given by the Bell states up to total phases). We no- 
tice that the Aj are constrained by Tr a = 0, then they 
sum up to zero. It is always possible to rearrange the 



coefficients as 



^ 9z ^ 0, using their 
properties of symmetry {6j f ^ ^j) £^nd periodicity 
{9j Oj + kj^, kj € Z) . The relation between the two 
families of parameters is given by 



A3 = —Or. — 



X <Jy 



A2 — + Oy 

^4 — Ox + Oy — 



(4) 



The aforementioned entanglement capability of the inter- 
action is defined as h = 0^ + Oy (e.g. see [5|). 

Notice that the decomposition ([2]) justifies the expres- 
sion ([1]) ; the non-local part is fully parameterized by three 
real constants, independent of the local actions. 

Dynamical evolution. — The most general Hamiltonian 
terms are given by 



Hi ^ ujifi ■ a 
Hi 



H2 ~ LU2I ®'m-a, (5) 

CyO-j (g) (Tj, 



where i and j range over {x,y,z}, oji and 0^2 are the 
characteristic energies of the two subsystems, n and m 
are real unit vectors, a is the vector of Pauli matrices and 



the coefficients 



form a real matrix. Without loss of 



generality, we will consider representations of the Pauli 
matrices, for the two subsystems, such that this matrix 
is diagonal, Cij = CiSij. 

The local actions consist of arbitrary preparations of 
n, m, bJi, and UJ2- In order to have manageable expres- 
sions, we assume that only the characteristic energies can 



be modified, and we fix n = m = (0,0, 1). We further 
consider ^ 0, i = 1,2. The Cartan decomposition 
of the unitary transformation Ut acting on the system is 
written as 



-iHrt _ 



LtAtKt, 



(6) 



and the dependence on time is made apparent. We are 
interested in the relevant contributions for the entangle- 
ment generation, that is At and Kt- It is possible to 
compute 



Ut = 











-icz t 



e"'''ip2(t) e^^-'ifisit) 
-*'=-V4(t) e- 



^V4(i) 





where 

fi{t) = cosfiit 

Vz{t) = 
with frequencies 



(7) 



I I — ^ J smiZii, 

cosli2r — * ( — t: I smsi2C, 

"2 ^ 



(8) 



^1 = \l [Cx - CyY + {lOi + UJ2Y , 



no 



\j [Cx + + (Wl - ^2)^ 



(9) 



In order to find the terms of the decomposition ([6]), it 
is convenient to represent all the operators in the magic 
basis, in which the local contributions become orthogonal 
matrices Lt and Kt, and the non-local part At is diagonal. 



Ut = LtAtKt 



(10) 



Since [// Ut — Kt AfKt, it is possible to determine At 
and Kt by considering the eigenvalues and eigenvectors 
of this operator. Applying the same procedure to UtUf 
it is possible to find also Lt, however this contribution is 
not relevant for the purposes of this paper. Denoting the 
eigensystem as 



U[Utv,{t)=eUt)v,{t), i = l,...,4, 
we obtain 



(11) 



£3.4V 



(12) 



where we have defined two real-valued functions, 



Xi{t) = l-2(^^j sm'n2t, 
X2{t) - l-2(^^)%in2l]it. 



(13) 
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Since —1 ^ 1 (« = 1,2), it is possible to check 

that \£i{t)\ = 1, i = 1,...4. The corresponding (non- 
normahzed) eigenvectors are given by 



V3A{t) = (Ai2WTJ^2(i),0,0,lf , 



(14) 



where 



/ii(t) = — — — cotf22t, /^2(i) = — — — cotfiit, 



OJl — UJ2 



(15) 

and Vi{t) = + y.f{t), for i — 1,2. The normahzed 
eigenvectors form the rows of the matrix Kt- 
The non-local part in the magic basis is given by 

At = diag{ei{t),i = l,..A), (16) 
and, following equation ([3]), it is possible to write 

£,(i) = e-'^^(*), j = l,...4. (17) 
Therefore the Cartan coefficients are given by 



Ai,2(i) = C2t± -arccosx2(i), 

^3,4(0 = -c^t T ^arccosxi(07 (18) 



and the expressions for the 0i(i), i = x,y, z, can be ob- 
tained inverting Q. These parameters characterize at 
every time the non-local contribution to the dynamics. 
Unlike the arbitrarily fast LU control case, they are time- 
dependent real functions that contain an explicit depen- 
dence on the local controls through the parameters fli 
and ri2, unless [H1+H2, Hj] = 0. This condition is satis- 
fied only in two cases: either uji — uj2 — 0, or Cx — Cy — 0. 
In both cases Aj (t) are linear in t and independent of the 
local controls, since the local dynamics and the interac- 
tion become independent processes. 

Finally, the relevant contributions can be rewritten in 
the original basis. The local term can be cast in the form 







ojj)2(t) 






-ir)i(t) 





oir;i(t) 






-im(t) 






(19) 



where 



Vi {t) - (t) ' 

The non-local contribution has the form 



i = 1,2. 



(20) 



At 



£2{t)+e3{t) 



£2(0 -£3(0 



£2it)-e3{t) 

ei{t)+ei{t) ei{t)-ei{t) 
£i(t) -£4^ ei{t)+e4{t) 



(21) 











£2(0 +£3(0 



Entanglement capability and optimal input states. — 
We are now able to derive the entanglement capability 
h{t) of Ut and study its properties. Considering p8)) . it 
is possible to obtain 

Ht) = I Kt) + ^ s{t) arccosx^(t)(^), (22) 

where k(t), s(t) and i{t) are piecewise continuous func- 
tions with values in Z, { — 1,1}, and {1,2} respectively, 
such that the hierarchy relations among the 9j{t) are 
fulfilled at every time t. The specific form of these 
functions is not relevant for the purposes of this work. 
The entanglement capability is a continuous function 
in t, with ^ h{t) < f , and it depends on the lo- 
cal parameters uji and UJ2 through the functions Xi(^), 
i — 1,2. The extremal points for Xi{t) can be found 
from dtXiit) = d^^Xi{i) = d^^Xi{t) = 0, whose relevant 
solutions are 



where the upper sign holds for i = 1, the lower for i = 2, 
and fc g Z. A typical dependence of h{t) on lji and lj2 for 
a critical value of t is represented in Fig. [1] The compli- 
cate peak-valley pattern is determined by the extremal 
points of the functions Xi{^)^ according to ([^^ . From 
it can be seen that the initial preparation of the 
local parameters has a strong impact on the evolution of 
the entanglement capability. In general, a sudden change 
is observed when uji and L02 match. This is consistent 
with the behavior of the purification process discussed in 
[l^l. In fact, the entangling capability of the evolution is 
fundamental in indirect control schemes. 

An initial state |'0)o is transformed by the evolution 
in the usually entangled state \'4')t — Ut\'ip)o- The max- 
imal attainable entanglement depends on the entangle- 
ment capability as well as on the initial state. It is possi- 
ble to characterize the set of the optimal input states by 
solving the equation 

\i^)^ = KtA\\^), (24) 



FIG. 1: Dependence on uJi and 0^2 of the entanglement capability 
h(t), for the Ising interaction Cx = 1, Cy = = 0, and t = ^n. 
The main peaks of h{t) correspond to uii = aJ2. The other peaks 
are related to the extremal points of Xi(*) with respect to loi and 
UJ2, and t fixed. 



FIG. 2: Dependence on uji and 0^2 of the entanglement capability 
h{t), with Cx = —0.2, Cy = 0.5, Cz = —0.8, t = 3.7, and n and 
m are the unit vectors corresponding to (1,-3,2) and (—2,1,4) 
respectively. The condition u!i = L1J2 still plays a special role. 



where |^) is an arbitrary maximally entangled state, Kt 
and At are expressed in ((T^ and (pij) respectively, and 

Conclusions. — When a two-qubits system is manipu- 
lated via local controls, it is usually assumed that these 
actions can be performed in an arbitrarily small time. 
Under this hypothesis, there is a clear separation between 
local and non-local contributions in the dynamics. The 
entangling capability of the evolution is a constant em- 
bodying the interaction content of the dynamics, and ev- 
ery initial factor state can produce the maximal amount 
of entanglement, since it can be instantaneously trans- 
formed, by fast local actions, in the optimal input state 
for the entanglement generation. 

In this paper we have explored the entanglement gen- 
eration in the opposite situation, where the local actions 
are fixed during the evolution. We have found that, in 
this regime, the local and non-local parts of the dynam- 
ics are strictly interconnected. Considering a particular 
control model, we have found the expression of the en- 
tanglement capability on the evolution time t, and on 
the local parameters ui, L02. In particular, h[t) reaches 
its global maxima periodically in i, under the condition 
LOi =0^2. Only some selected input states maximize the 
entanglement production. 

If less restrictive control models are adopted, in general 
it is not possible to obtain simple analytical expressions 
for the entanglement capability. However, there is strong 
evidence that the main features of the behavior of h{t), 
described in this work, do not depend on the particular 
choice of n and m. In fact, we have numerically com- 
puted h{t) for a large sample of Hamiltonian operators, 
with randomly distributed values of c^;, Cy, n and rfi. 
We have always found results that are consistent with 
the analysis presented in this work, in particular the in- 



surgence of the peak of h{t) when the two characteristic 
energies match. An prototypical example of these simu- 
lations is presented in Fig. [21 that exhibits the behavior 
of h{t) for a particular evolution with n ^ rfi. The in- 
crease of h[t) in correspondence of wi = lj2 is not an 
artifact of the choice of the Hamiltonian terms assumed 
in this work, it is rather a general phenomenon in the 
two-qubits system. 
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